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Abstract 

We consider the one-dimensional totally asymmetric simple exclusion model (TASEP model) with open bound- 
ary conditions and present the analytical computations leading to the exact formula for distance clearance dis- 
tribution, i.e. probability density for a clear distance between subsequent particles of the model. The general 
' relation is rapidly simplified for middle part of the one-dimensional lattice using the large N approximation. Both 

the analytical formulas and their approximations are successfully compared with the numerical representation of 
the TASEP model. Furthermore, we introduce the pertinent estimation for so-called spectral rigidity of the model. 
\^ ' The results obtained are sequentially discussed within the scope of vehicular traffic theory. 
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PACS numbers: 02.50.-r, 05.45.-a, 89.40.-a 



1 Introduction and motivation 

The effect of queuing is one of commonly appearing phenomenons in nature. It can be disclosed in microbiological 
systems, societies of animals, computer networks, public transport systems, parking cars manoeuvres and many oth- 
ers. In the recent years, queuing has attracted attention of many physicists and mathematicians. It is well known 
f*^; ■ that many of above-mentioned dynamics systems belong to the same class of mathematical tasks. A convenient way 
• | how to investigate such systems in detail can be found in space-discrete and time-continuous one-dimensional models 
■ based on asymmetric exclusion processes. Indeed, the family of asymmetric exclusion processes is as wide-ranging 
| as family of its applications. Concretely, the ASEP-models have been successfully used for description of protein 
t-H . synthesis [I], [2], polymers in random media [3], fluctuations in shock fronts [4], [5], [6], [7], gel electronic [8], in 
; molecular biology [9] and finally in physics of traffic [TDJ, [11], [12], [13] or [14]. 

In this article we are focused predominantly on those aspects of cellular modeling having a vehicular-traffic 
jj^ • interpretation. With respect to the fact that macroscopical behavior of TASEP (represented for example by the 
fundamental relation between the flux and density of the model) corresponds to that detected in the freeway samples 
(compare Figures 1 and 4 for illustration) it is meaningful to consider the TASEP-model as an appropriate choice for 
one-lane traffic simulator. However, besides the confrontation of macroscopic quantities it is indispensable to make 
a comparison of microscopic quantities. Recently, the extensive investigations of traffic data ([18], [19] , [20J, [21J, 
[22J, |23J, or [24j) provide a good insight into the microstructure of traffic samples. It has been demonstrated in 
some of the above-mentioned articles that inter-vehicle gap statistics (clearance distribution) in real-road traffic can 
be very well estimated by the one-parametric family of functions 

p(r)^e(r)e^e- Br , (1) 

where 

5 = ,+ ^— , (2) 

(3) 

We remark that the functions 6(.t), IC\(x) represent the Heaviside's step-function 
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Figure 1: The Empirical Flow-Density Relation and Relevant Schematic Represen- 
tation. The dependence of traffic Flux J on traffic Density g for realistic traffic 
flows uses to be compared to a mirror image of the Greek letter A (see |15j or [16J). 
The relation J = J(g) extracted from the induction-double-loop-detector data is 
visualized in the inset. 



and the modified Bessel's function of the second kind (Mac-Donald's function), respectively. Furthermore, the 
distribution p(r) fulfils two normalization conditions 

/ P(r)r = l (4) 

JR 

and 

/ rp(r) r = 1. (5) 
Jr 

The latter represents a scaling to a mean clearance equal to one. We add that the one and only parameter v is related 
to the traffic density g (see [25j for details). Roughly speaking, such a parameter (called speculatively as mental strain 
coefficient) reflects a rate of psychological pressure under which the car-drivers are during driving manoeuvres. The 
changes of traffic status from free flows to congested flows and vice versa are accompanied by the adequate changes 
of mental strain coefficient v 1 i.e. by the adequate changes of clearance distribution. The chosen representatives of 
the relevant analysis are displayed in the Figure 2. Here one can detect the basic probabilistic trends of distances 
among succeeding cars. Whereas in the region of small densities (free traffic regime) the relevant probability density 
is exponential essentially (which fully corresponds to the fact that cars interactions are negligible) the distribution 
p(r) is rapidly changing if congested data are observed. In this case the stronger mutual interactions among vehicles 
lead to the hardcore repulsions in the system, which results in fact that 

lim p(r) = 0. 

r->-0 + 

As discussed in [25j the intermediate region of metastable traffic states shows a substantial growth of parameter v. 
This is influenced by the fact that driver, moving quite fast in relatively dense traffic flow, is under a considerable 
psychological pressure. After the transition from free to congested regime the pressure momentarily declines because 
of decrease in mean velocity. Finally, if the traffic flow becomes denser and denser the mental strain coefficient v is 
increasing further. This finally culminates by the creation of stop-and-go traffic waves. 

The maim goal of this article is to describe some correspondences between the microstructure of asymmetric 
simple exclusion model and real-road traffic. For this purpose we will compare (in the first part of this work) the 
relevant clearance distributions, i.e. probability densities for clear gap among all pairs of succeeding particles/vehicles. 
Subsequently we will analyze the associated spectral rigidities of both systems and will discuss their similarities. 
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p ~ veh/km 




Inter-particle Gap r 

Figure 2: Clearance Distribution for Real-Road Traffic. The curves visualize the 
analytical estimations ([I]) of the inter-particle statistics obtained for the densities 
indicated in the legend. 



2 The totally asymmetric simple exclusion process (TASEP) with 
open boundaries 

Consider a chain of length N containing TV equivalent cells and define three fixed parameters a,/3,p E [0, 1]. Let each 
cell to be either occupied by one particle or empty. During the infinitesimal time interval dt each particle hops to the 
immediate site (in the defined direction) with probability pdt if the target site is empty. In the opposite case (if the tar- 
get site is occupied) the particle does not change its location. If the first cell (£ = 1) is empty a new particle is injected 
into the chain with probability adt. Similarly, if the last cell (£ = N) is occupied the relevant particle leaves the chain 
with probability (5dt. This definition can be generalized if needed, however, for all purposes the original formulation 
is fully sufficient. Furthermore, we will use (without loss of generality) the re-scaled variant of the model where p = 1. 



adt 




pdt 



ooo»oo»oo«o 




Figure 3: Definition of the TASEP. The configuration visualized can be identified 
with the vector C = (0, 3, 1, 2, 1, 2, 1, 1, 2). 



The sketched rules define a simple driven lattice model of one-lane traffic whose elements are hard-core-repulsed 
by "implicit forces" derived from the above-mentioned exclusion rules. A great merit of such a model lies in the 
exact solvability of the associated steady state. Indeed, two alternative methods for exact solution of TASEP are 
outlined in [26j, [27]. For intentions of this research we consistently use the methods based on the matrix formulation 
introduced in [27]. According to this matrix method we consider the matrices D, E (infinity-dimensional, in general) 
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and vectors (w\ and \v) satisfying the following algebraic rules 

DE = D + E, (6) 

{w\E=hw\, (7) 
a 

D\v) = (8) 
Although each configuration of the model can naturally be described by the binary sequence 

T=(t 1 ,t 2 ,...,t n )g{0,1} N , 

it is profitable to use the following convention. Let the symbol C = (nx,mx, n 2 , m.2, • ■ • , n q , m q ) denote the configu- 
ration of the TASEP-chain such that (starting from the initial cell of the chain) n\ is the maximal number of empty 
cells (i.e. the (ni + l)th cell is occupied), mi is the maximal number of occupied cells (starting from the (m + l)th 
cell), when the (m + mi + l)th cell is empty, and so on. Thus, J2i n i + J2j m j = N. The previous definition 
(visualized for lucidity in the Figure 3) corresponds unambiguously to the relevant arrangement of the particles inside 
the TASEP-chain, which means that there exists some bijection between T and C. Moreover, it has been proved in 
|28j that the steady-state probability of an arbitrary configuration C reads as 

1 (w\D ni E mi D n *E m2 ...\v) 

r ss (C) 



Z N (w\v) 
where the constant Z^ assures the proper normalization, i.e. 

]Tp ss (c fc ) = i 

c k 

for k running over all possible configurations. In fact, such a probability can be calculated directly using the relations 
ecJU j?}, and ([8]), i.e. without knowingthe specific representation of matrices D, E and vectors (w\ and \v). Moreover, 
let the symbol 

, _ J 1, i — th site is occupied 
^ ' [ 0, i — th site is empty, 

represent (for the fixed configuration T) the binary functional ojj- : N n- {0, 1}. Let the symbol 

= {n 1 ,mi,£ 1 ,n 2 ,m 2 ,i2, ■ ■ • , n q ,m q ,£ q ) 

denote the set of all configurations such that ni is the maximal number of empty cells (i.e. the (rii + l)th cell is 
occupied), nil is the maximal number of occupied cells (starting from the (ni + l)th cell), when the {ni + mi + l)th 
cell is empty, and £i is the number of cells (behind the (?ii + mi)th cell) which can be arbitrarily occupied or not. 
Precisely speaking, the symbol 2£ corresponds to the set 

{(wr(l), wr(12), ■ • • , <*r(N)) e {0, 1} N : w r (l) = u T {2) = ... = wr(ni) = A 
wr( n i + 1) = u T {ni + 2) = . . . = bJr{ni + mi) = 1 A 

wr( n i + m i + £i + 1) = wr(fii + mi + 1\ + 2) = . . . = u-j-(ni + mi + £i + n 2 ) = A 

UJri n i + mi + £i + n 2 + 1) = wr( n i + m i + + n 2 + 2) = . . . = wr( re i + mi + £i + n 2 + m 2 ) = 1 A 

Wr( n l + mi + £i + ri 2 + m 2 + £ 2 + l) = . . . = uJr{ n i +mi+ £i+ n 2 + m 2 + £ 2 + n 3 ) = A . . . j. 

As follows from the results published for steady state of TASEP-model (see for example |28J or |29j ) the probability 
for finding the system in one of configurations included in 5£ is 

1 (w\D ni E m i{D + E) ll D n *E m *(D + Ey-> ...\v) 
Z N \w\v) 

As a direct consequence of this assertion we easily deduce that if X = (0,0, iV) then 

P ss(t5 n = 1 (w\D°E°(D + E) N \v) = 1 (w\(D + E) N \v) 
Z N (w\v) Z N (w\v) 
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As 3£ = (0,0, N) includes all permissible configurations of the system (and therefore P SS (J£") = 1) one can trivially 
calculate the value of the respective partition sum as 

(w\(D + E) N \v) 
<&n = — 



(w\v) 

Using the lemma [5] and notation (f33|) we can assert that 



1 N /our |ii ™ N rn 

z - = t-tt E e Ew ^» = E B -.» E ■ 



(lulu) ^ AT!(iV-m)! ^ x 1 ^ ' ^ of 8 

x 1 ' m=l v ; i=0 m=l i=0 ' 

Hence 



Lm=l B JV,m ^ 8-i - a -i « 7* P 



Z)m=l B »,™ mQ ™ 01= B. 



(10) 



We conclude that by means of the formulas ([9]) and (fT0|) one can enumerate the probability of an arbitrary 
steady-state-configuration of TASEP. 

3 One-dimensional representation of associated matrix algebra 

As demonstrated above, the steady-state probability distribution of an arbitrary arrangement of TASEP particles is 
derived from the matrix algebra ([6]) - ([8|). Such a quadratic algebra is formulated for two matrices D and E which 
are associated to the particles and holes respectively. As proven in [27] the matrices fulfilling the rules {[G]) — dSD are 
non-commuting and infinity-dimensional, in general. However, for special choice of parameter a and B the matrices 
E and D can be commuting, i.e. DE = ED. Under this condition the equation ([6]) leads to the equations 

D + E = DE = ED, (11) 

jj\v) + E\v) = Ejj\v) 

~{w\v) + -(w\v) = (w\^-\v) 
p a ap 

a + 8 = l. (12) 

Thus, if the condition (j!2p is guaranteed the commuting matrices D and E are one-dimensional, i.e. they are 

represented in fact by the numbers E = aT 1 and D = B^ 1 . Consecutively, we can choose (w\ = \v) = 1. Above 
that, as the partition sum is now reduced to 

(w\(D + E) N \v) fl^ N 

£>N = - 



(w\v) \ct8 
the steady-state probability (for the set of configurations reads as 

P SS (,T) = (a/?) \j) [^) =a 8 = a (1 - a) . 

4 Macroscopic characteristics in mean-field approximation 

For intentions of this article we briefly summarize the known results on macroscopical behavior of TASEP-model. We 
are concentrated predominantly on the bulk density and flux in mean-field approximation. Concretely, the average 
density of particles occurring inside the ith cell can be depicted as a probability for finding the system in the certain 
configuration chosen from 3£ ei = (0, 0, i — 1, 1, 0, N — i). That means 

W _ p (<r {w\{D + Ey- l D(D + E) N ^\v) 
g t -f ss [^ gi )- ( w \(D + E) N \v) 

Similarly, the flux through the ith site is defined as 

tW _ p (»■ \- (w\(D + Ey^DE(D + E) N -^\v) _ Z^v_i 
" ss[ Ji) (w\(D + E) N \v) Z N ' 
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where S£j i = (0, 0, i — 1, 1, 1, N — i — 1). Using the large N approximations 

N+l 



@a(l-2a) 


( 1 


P—a 




/3a(l-2/3) 




a-P 




Pa 


1 


y^(a-p) 


(2a-l) 2 


a 2 


" 4 N + l 


J^(2a-1) 3 N 3 / 2 


(l-2<*) 2 


N 


(1-a) 2 a N (l-a) N 


2P 


4™ 


y^F(2/3-l) N 1 / 2 



N+l 



7V-V2 



a < | A a < /3 

/9 < | A P < a 

a > \ A P > \ 

a = P>\ 

oi = p <\ 

a. = \ < P 

OL~P=\ 



(13) 



derived in |27j we ascertain that (for an arbitrary cell being far from the boundaries of the system) the relevant 
bulk-density is given by 

\ A P^\, 

a a < i A P > a, 

1 - p P < ± A P < a, 
a a + P = 1. 

Regarding the TASEP flux J it is trivial to show that it is independent on i and 



g(a,p) 



(14) 



J(a,p) = 



a(l — a) 
P(l-P) 
a(l — a) 



a ^ | A P > i, 
a < | A > a, 
/3 < | A /3 < a, 
a + p = 1. 



(15) 



After comparing the results (fl3|) and (|15p we close this section by the assertion that the fundamental diagram of 
TASEP-model (see Figure 4) is described by the equation 



J(q) = g(l - g). 



(16) 




Figure 4: The Flow-Density Relation for TASEP-model. The curves display the 
function J = J(g) calculated for ordered-sequential update of ASEP-model with 
periodic boundary conditions (for details on TASEP update-procedures, please, see 
|17|). Parameter p represents the rate for particle transition from the actual cell 
to the following one. The bold curve displays the theoretical relation (j!6p valid for 
time-continuous TASEP. 



5 Distance clearance distribution - a special case 

In the main part of this article we are aiming to derive an exact analytical formula for the so-called distance clearance 
of TASEP-model, i.e. probability density for clear distance (measured in cells) among the subsequent particles of the 
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model. Such a probability density can be (in general) depending on the position i of a monitored particle. Hence, 
the relevant distribution pi(k) is depending on i as well. The unnormalized probability of a fccell-gap behind the ith 
occupied cell is denoted as pi(k) = P ss (0, 0, i — 1,1, k — 1, 0, 1, 0, N — i — k) and can be calculated via 

,,, 1 WD + B r>DEt->DiD + E )-»M p<i<w _ 0i (17) 

Z N (w\v) 

because the first i— 1 cells could be randomly occupied or vacant, the zth cell should be occupied, k — 1 succeeding 
cells should be vacant as well, the (i + fc)th cell should be occupied, and finally the tail of the chain could be 
occupied/vacant at random. The probability density for clearance is then p[ N \k) = J\fpi N \k), where the constant 
M assures the proper normalization computable via the condition 

N-i 



k=l 



If the associated matric algebra is one-dimensional, i.e. under the condition (fl2|) , then (|17jl changes to 

, 1 x N-2 , 1 x N+l-k 

pf\k)=M - hi) (l<k^N-i). (19) 



a) \/3 

Precisely speaking, the derivation of the last formula using the relation (fT7|) is legitimate for k 3, since for k = 1, 2 
the relation (j!7p has a slightly modified form. However, it can easily be verified that relation (|19p holds true for all 
positive k. The normalization pre-factor Af can be quantified either from condition (fl8|) or from the equation 



x _ (w\(D + E) l ~ 1 D(D + E) N ~ l \v) - (w\(D + E) l ~ l DE r 



that takes into account all possible configurations with the occupied ith cell and rejects those configurations having 
the empty tail (behind the zth cell), when all the cells i + 1, i + 2, . . . , N are vacant. This leads to 



JV-l / 1 \ JV 

TV -1 = f - 1 



which yields (if (fl9|) is applied) the result desired 

The large N approximation adjusts such a expression to 

Pi (k)= lim pf\k) =a(3 k -\ 

N— f+oo 

Note that this result is independent of the cell behind which the clearance is measured. As a consequence, we obtain 
the formula 

(*) = p k Pi (k) = a p k^ = = i 

for mean clearance of TASEP-model. If reformulated the clearance distribution reads as 

p(k) = g(l - o) k -\ 

where the relation g = a is adopted. We remind that these outcomes are valid under the condition a + /3 = 1 only. 

6 Distance clearance distribution - a general case 

To derive the exact formula o 
product in (fT7|) into the form 



To derive the exact formula of p\ N \k) for an arbitrary choice of parameters a, ft we will first rearrange the matrix 



(D + E) i ~ 1 DE k ~ 1 D(D + E) N - l ~ k = J2 a m>n E m D n . 
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Then it will be hold 

(w\(D + E) i ~ 1 DE k ~ l D{D + E) N - l - k \v) = ^ a m , n cT m $- n . 
Using the notation ([33]), lemma[4] and lemma[5]we get (for k ^ N — i) 

N-i-k m N-i-k m / £ \ 

D(D + E) N -*- k = ]T B N + k , m y £ l DE t D m - t = £ B w _,_ fe , m £ £>»-<+i + £ S*D"»-' 

m=l £=0 m=l f=0 V 2=1 > 

Furthermore, by means of lemma[5]we obtain 

N-i-k m / £ \ 

DE k ~ 1 D(D + E) N - l - k = J2 B N-i-k, m DE k - x D m - g+1 + DE z+k ~ 1 D ri 



•\m—£ 



m=l 1=0 \ 2=1 / 

N-i-k rn / fc-1 £ 2+fe-l 



J2 B N _ t _ k , m Y,[D m - e+2 + j2 EWDm ~ e+1 + iDm ~ e+1 +J2 E EWr>r ' 



m=l £=0 \ w=l 2=1 to=l ) 

For i ^ 2 we apply the lemma[4]in the product (D + E)*~ . That provides 

-t-1 p N-i-k 



(D + Ef-'DE^DiD + E) N -*- k = EE E E fe, m B 4 - 1 . p (^^-" +m - £+2 + 

p=l g=0 m=l £=0 

fc-1 £ 2 + fc-l 

_|_ £jq jjp-q-^iu jjm-l+1 _|_ ^g^P _ 9+ m_£ + 1 + D p ~ q E w D m ~ e 

U) = l 2=1 10 = 1 

From the lemma[6]we derive that (for q < p) 

DP - qEW = y (p-g-a + w- 1)1 P yi (p- q -b + w -l)\ = ^ ™ 
^ (p-g-l)I(ttf-o)! ^ (p_ g _6)!(«,_i)! £j 

where we use, for convenience, the notation 

(p — q — a + w — 1)! (p — q — b + w — 1)1 

Al = (p - q - l)\(w - a)\ ' A2 = (p-g-fc)!(iu- 1)!' 

This leads to the final formula 

(D + E) i ~ 1 DE k ~ 1 D{D + E) N - i ~ k = 

i—l N—i—k m p— 1 / k— 1 ID 

= E E E E BiV-i-fe,mBi-l,p £9DP-9+"'^+ 2 + £ E Q D P-9+m-e+l + ^ g Ai^ <?+a J D m - £+1 + 
p=l m= l £ = o <j=0 \ io=l o=l 

k-lp-q £ 2+fe-l uj £ 2+fc-l p-g \ 

+ EE A2 ™ m ^ +1+b +E E E A i^ 9+a ^ m ^ + E E E A2 ™ m ~* +6 + 

u) = l a=l 2=1 to=l a=l 2 = 1 io=l 6=1 / 

i-1 N-i-k rn / fc-1 £ z+fc-1 N 

+ E E E BjV-i-fc,mBi-l,p E p D m ~ e+2 + E w+v Dm -£+l + eE p D m-t+l + ^ ^ E w+p Dm -£ 

p=l m=l £=0 V io=l 2 = 1 io=l > 

Multiplying such an expression by vectors (u>| and |i>) gives the following dependence on parameters a, (3 

9 =0 v 7 u)=l ^ \a=l v 7 6=1 
z=l q=0 x 7 to=l \a=l v 7 6=1 



where we consider i 2 and k < N — i. The formulas for the remaining cases have been obtained analogically. 

. N—i — k m / -i \ m ~ ^ f i n k—X / 1 \ w 1 ^ z+fe— 1 ✓ - 



9 =0 v 7 io=l \a=l v 7 6=1 



^-«-G) ,+ s(=)'(J 



u;=l 

Now we are focused on the asymptotic clearance distribution p(k). As we are interested in the probability of a 
gap of the length k between two succeeding particles in the bulk, i.e. far from the boundaries, it is beneficial to fix 
the position of cells i, i + 1, , . , , i + k. Based on this strategy we denote the number of arbitrary occupied cells before 
the gap as m = i — 1 and the number of arbitrary occupied cells behind the gap by n = N — i — k. The limit N — > oo 
is then meant in the sense m,n — > oo. This approach assures that the gap is located far from the boundaries with 
respect to the lattice size N, which is necessary for the calculations below. 

Let us now investigate the behavior of the quantity p m ,n{k) = p^+i{^)< wnere N = m + n + k + 1. In the term 
of matrix-product-ansatz description (MPA) we focuss on 

~ 1 (wKD + ErDE^DjD + Eriv) 

Pm,n{ k ) = ~ ]—T\ ( 21 ) 

Z m +n+k+l (W\V) 

in the limit m,n — > oo. Firstly we derive the asymptotic form p(k) = lim„ l: „_ i . 00 p m . n {k) for k = 1,2 and then we 
get a recursive formula for p[k + 1), k ^ 3. Rearranging the matrix product in (|2T|) using (fTT| ) for k = 1 we get 

(D + E) m DD{D + E) n = {D + E) m D{D + E) n+1 - (D + E) m DE(D + E) n = 

= (D + E) m D{D + E) n+l -(D + E) m+n+l . 

And hence, using the asymptotic form of fundamental dependence J(g) = g(l — g), it holds 

p m ,n(l) = & m+2) - J [n+m+2) => P{1) = Q-Q{l-Q) = Q 2 . (22) 
For k = 2 we analogically obtain 

(D + E) m DED(D + E) n = (D + E) m+1 D(D + E) n , 

and therefore 

p m ,„(2) = ^ + + ! n+2) j(™+«+ 3 ) p(2) = gg{l-g) = g 2 {l-g). (23) 

Moreover, for k ^ 3 it holds 

(D + E) m DE k - 1 D(D + E) n = (D + E) m DE k ~~ 2 D(D + E) n + (D + E)" 1 ^ 1 DE k D(D + E) n 

That leads to the formula 

p m -x, n (k + 1) = p m ,„(fc) - J m+ " +fc+1 p m ,„(fc - 1). (24) 
For an arbitrary fixed k £ N it holds that 

lim p m .n(k) = lim p OT _i n (fc) =: p(k). 

m : n— foo m,n— ^oo 

The existence of the limit p(k) for all k € N can be elementarily confirmed using mathematical induction by means 
of relations (|22p , (|23p, (|24p. Now we can apply the limit m,n->ooon the relation (|24p obtaining 

p(fc + i) = p(A)-e(i-e)p(ft-i), p(i) = ^ 2 , p(2) = g 2 (i - £). 
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It is easy to verify that the solution of this recursion is the relation 

P(k) = g 2 (l - g)*" 1 . 

The normalization constant Af is then 

and hence the clearance distribution p(k) = Afp(k) fulfils the relationship 

p(k) = g (i- e ) k - 1 . (25) 

All the time we kept in mind that the bulk density g is a function of parameters a, j3 given by (|14p. Thus, the 
dependence of the clearance distribution on these parameters reads 

a(l - a) fc_1 a < i A /? > a, 
(1 — /3)/3 fc " 1 f3<±A(3<a, 
a^- 1 a + 13=1. 



p(k;a,(3) 



7 Mean-field spectral rigidity of TASEP-model 

Description of traffic microstructure by means of inter-vehicle gap distribution p(r) is, as discussed above, usual in 
traffic theory. However, clearance distribution depicts the distance gaps between two successive vehicles only. Aiming 
to investigate the middle-ranged interactions among the cars it is necessary to find a mathematical quantity suitable 
for quantifying the level of synchronization for larger clusters of particles. This desired quantity can be found in 
Random Matrix Theory (see [30]) where it provides an insight into the structure of eigenvalues of random matrix 
ensembles. It is called a spectral rigidity. If reformulated within the bounds of traffic theory the spectral rigidity has 
the following interpretation. 



+ 




Scaled Inter-vehicle Gap 



Figure 5: Spectral Rigidity of Vehicular Traffic. Gray curves correspond to the 
rigidity A(L) analyzed separately in 85 density regions. The chosen results of 
relevant statistical analysis are picked out. Concretely, plus signs, squares, triangles, 
stars, circles, and diamonds represent the rigidities obtained for traffic data from 
the following density regions: [2,3); [6,7); [16,17); [26,27); [40,41); and [85,86) 
vehicles/km/lane respectively. For details, please, see the article [25] , 



Consider a set {r.; : i = 1 . . . Q} of netto gaps between each pair of the succeeding cars moving on an one-lane- 
freeway. We suppose that the mean gap taken over the complete set is re-scaled to one, i.e. 

Q 

5>i = Q. (26) 
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Dividing the interval [0, Q) into subintervals Uk — l)L, fci) of a length L and denoting by rife(L) the number of cars 
in the fc— th subinterval, the average value n(L) taken over all possible subintervals is 

1 IQ/H 

" (l) = Iq7ZJ £ »*( L ) = L > 



fe=i 



where the integer part [Q/L\ stands for the number of all subintervals Uk — 1)L, kL) included in the entire interval 
[0,Q). We suppose, for convenience, that Q/L is integer, i.e. [Q/L\ = Q/L. The spectral rigidity A(L) is then 
defined as 

T Q/l 

A(L) = -J2(n k (L)-L) 2 

w fe=i 

and represents the statistical variance of the number of vehicles moving at the same time inside a fixed part of the 
road of a length L. For the one-parametric family of distributions ([TJ there has been proved in the article [25] that 
the associated spectral rigidity has a form 

A(L) = X L + 1 + 0(L- 1 ), (27) 



where (together with the relation ([2])) 



and 



x = xiy) = ; 28 



6VBl} + Bv(2l + ABv + lGVBiA 

7 = 7M = =4 29 

24(1 + y/Bv) 

It means that the rigidity is a linear function whose slope \ is depending on the mental strain coefficient v (briefly 
described in the section [l]). Approximately, the changes of % (with respect to the traffic density) can be described 
as follows. For small densities the slope \ in relation A(L) w \L + 7 is close to one as expected for statistically 
independent events. Nevertheless, if the traffic density increases the interactions among vehicles strengthen, which 
results in a descent of the slope \- The more detailed insight into the realistic behavior of A(L) is demonstrated in 
the Figure 5, where the experimental results of "traffic spectral analysis" are plotted. 

At the moment our immediate goal is to derive the analytical formula for the spectral rigidity A of totally 
asymmetric simple exclusion model. Firstly, note that the mean clearance calculated for the distribution (|25]l is 



00 

{k)=Y J kQ{l-Q) k - 1 = 



fc=l 

which means that in the end of our computations there will be necessity to revise the results in accord with the 
general definition (|2bp. Now, denote by ni{s) the probability that there is exactly s particles inside the fixed I— cells 
region of TASEP-chain. As understandable, such a probability reads 

n/(a)= (a = 0,1,...,*). 



Since 
and 

(s 2 ) 



(^E(^r^^-^ = ^ 



the corresponding statistical variance is of a form 

a 2 = (s 2 ) - (s) 2 = q£(1 - g), 
which leads (after the transition to the re-scaled headways L = g£) to the final formula for the spectral rigidity 

A(L) = (1 - q)L. (30) 

Obviously, the slope \ = 1 — Q of the rigidity A(L) decreases linearly with increasing density of the TASEP-particles, 
as confirmed by the numerous numerical tests visualized in the Figure 6. This result is very similar to the behavior 
revealed in realistic traffic samples. 
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Figure 6: Spectral Rigidity of TASEP-model. Signs represent the numerically- 
obtained spectral rigidity A(L) for TASEP with parameters a,f3 indicated in the 
legend. The curves display the corresponding rigidity calculated via the exactly- 
derived formula ([30]) . 



8 Summary and conclusions 

Using the MPA formalism we have derived the universal analytical formula for distance clearance among subsequent 
particles of TASEP with open boundary conditions. Such an exact formula has been substantially simplified for 
middle segment of sufficiently long TASEP-chains and consecutively compared to the empirically-obtained clearance 
distribution of vehicular streams. As uncovered, the interactions among the TASEP-particles are extremely short- 
ranged and hard-cored, which leads to the fact that the continuous (and re-scaled) alternative for associated clearance 
distribution (|25ji reads as 

p(r) = 9(r)e-'\ (31) 

Thus, in spite of the fact that the TASEP-particles are repulsed such interactions are significantly different from those 
detected in realistic traffic flow. Anyway, the traffic clearance is of the form (TjQ), which means that distribution (f3T|) 
corresponds to the traffic clearance distribution only if the traffic density is small {y ps 0). But on the other hand, 
if the spectral rigidities are compared one can detect a good agreement between the TASEP-rigidity (f30|) calculated 
analytically and the same quantity analyzed in traffic data (compare Figures 5 and 6). 

To conclude, although the traffic model based on the TASEP-rules shows similar macroscopical effects as systems 
of moving vehicles it is definitely not suitable for modeling the microstructure of traffic (especially in congested traffic 
regimes). 
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9 Appendix 

Lemma 1. Let m,n,a be the arbitrary natural numbers. Let a ^ m. Then 

E l (n — i + m— 1\ / n + m — a\ . . 

\ n-1 = \ n l { > 



Proof. We will proceed using the mathematical induction. Let m = a = 1. Then for an arbitrary neN 

^(n-i\ = (n-l\f 
,11-1/ V n — 1 / \n. 



i=l 
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n- 1 



□ 



(33) 



Now let the relation (|32|) hold for all a ^ m. Then for a < m + 1 

m+l 

/ n — i -\- m \ x ^ / n — 

tj — 1 / V 7i — l / ' V n — 1 / v n 

/ i—a x 7 x ' x 

The case a = m + 1 can be proven analogically to the case m = 1, a = 1. 
Lemma 2. If using the notation 

f p(2n- P -i)! < d < n 
g _ J n!(n-p)! ^ ^ ' 

[O otherwise, 
it can simply be verified that for n > and p ^ 

Lemma 3. For arbitrary square matrices D, E fulfilling DE = D + E and for n 

ra+1 

D n (D + E) = D n+1 + D n ~ 1 (D + E) = Y £>« + E. 

9=1 

Lemma 4. Assume that D, E to be arbitrary square matrices fulfilling the relation DE = D + E. Let n be the 
natural number. Then 



(D + E) n = Y B„, m E E9Dm ~ q - 

rn—1 q—0 



(34) 



Proof. Again, we will use the mathematical induction with respect to n. Let us assume that (|34p holds. Then with 
help of lemma[2]and lemma[3]we elementarily deduce that 

n m n m / m—q+1 \ 

(D + E) n+1 = Y B n,?nY, EqDm ~ q ( D + E )= E B »+l.™+lE [E q+1 +E1 J2 D " - 
m—1 q—0 m— 1 g— \ p— 1 / 

n m / m— g+1 \ n+1 m 

- E B «+i,m+2 E + S " E ^ = B «+2,2(^ + E f + E B„+i, m £ E«D m -i = 

rn—1 q—0 \ p—1 / m— 3 q—0 

n+1 m 

.7=0 



m— 1 



□ 



Lemma 5. Assume that D, _E to be arbitrary square matrices fulfilling the relation DE = D + E. Let n G N. Then 

n 

DE n = DE"- 1 + E n = D +Y EW - (35) 

Lemma 6. Assume that D, E to be arbitrary square matrices fulfilling the relation DE = D + E. Let m,n be the 
natural numbers. Then 



i=l v 7 j=l 



7i — j + m — 1 
rn—1 



(36) 



Proof. We will use the mathematical induction. For n = 1 and an arbitrary m £ N the expression (|36|) reduces to 
the expression (|35|) . Let us now assume that (f36|) holds true for some n E N. By means of lemma[T]and lemma[5]we 
get 



D(D n E m ) = Y 



i=i 



n — i + m — 1 



ri- 1 



71+1 



3 = 1 



u;= 1 
m m 



E 

3=1 



"-i + ^-A^+i = 



771 — 1 



-U;— 1 7—71! 



71 — 7 + 771 — 1 
77 - 1 



n+1 



3 = 1 



77 — j + 771 
777 — 1 



□ 
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Figure 7: Clearance Distribution in the Phase of Low Density I. We plot the inter- 
particle distance distribution calculated for TASEP with a = 0.1 and f3 = 0.2. 
The bullets visualize the clearance distribution p\ (k) enumerated via the exact 
formula ([20]) . while continuous curves visualize the large N approximations p(k) 
summarized by the relation (f25|). Above that, we also plot the relevant densities, 
which synoptically illustrates the reasons for discrepancies between (|20|) and (|25|). 
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Figure 8: Clearance Distribution in the Phase of Low Density II. We plot the inter- 
particle distance distribution calculated for TASEP with a = 0.3 and (3 = 0.9. 
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Figure 9: Clearance Distribution in the Phase of Maximum Current. We plot the 
inter-particle distance distribution calculated for TASEP with a = 0.8 and f3 = 0.9. 
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Figure 10: Clearance Distribution in the Phase of High Density I. We plot the inter- 
particle distance distribution calculated for TASEP with a = 0.9 and j3 = 0.3. 
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Figure 11: Clearance Distribution in the Phase of High Density II. We plot the 
inter-particle distance distribution calculated for TASEP with a = 0.2 and f3 = 0.1. 
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